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Abstract 

We consider the compute-and-forward paradigm with limited feedback. Without feedback, compute-and-forward 
is typically realized with lattice codes over the ring of integers, the ring of Gaussian integers, or the ring of 
Eisenstein integers, which are all principal ideal domains (PID). A novel scheme called adaptive compute-and- 
forward is proposed to exploit the limited feedback about the channel state by working with the best ring of 
imaginary quadratic integers. This is enabled by generalizing the famous Construction A from PID to other rings 
of imaginary quadratic integers which may not form PID and by showing such the construction can produce good 
lattices for coding in the sense of Poltyrev and for MSE quantization. Simulation results show that by adaptively 
choosing the best ring among the considered ones according to the limited feedback, the proposed adaptive compute- 
and-forward provides a better performance than that provided by the conventional compute-and-forward scheme 
which works over Gaussian or Eisenstein integers solely. 

Index Terms 

Compute-and-forward, physical-layer network coding, lattice codes, and algebraic integers. 


I. Introduction 

Compute-and-forward [|3 is a novel paradigm of information forwarding whieh allows relay nodes to 
compute and then forward functions of messages by exploiting the structure induced by the channel. The 
main enabler of the scheme in |l2l is the use of lattice codes from Construction A over Z the ring of 
integers. Later it was shown that compute-and-forward can also be performed over other rings of integers 
such as the Gaussian integers Z[i] and the Eisenstein integers iJyuj] [[3l-[|5l by constructing lattices via 
Construction A over the respective ring of integers. 

On one hand, since Z[i] and Z[a;] are instances of imaginary quadratic integers, it seems natural to 
extend the compute-and-forward framework to general rings of imaginary quadratic integers. On the other 
hand, since the role of the underlying ring of integers can be effectively thought of as being a quantizer of 
the channel and Z[a;] is already the best quantizer for C where the channel coefficients belong, it seems 
unnecessary to pursue compute-and-forward over other rings of integers. In this paper, we first seek to 
better understand the role of rings of algebraic integers in constructing good lattices. We then use an 
example, namely compute-and-forward with limited feedback, to demonstrate the benefits of performing 
compute-and-forward over rings of imaginary quadratic integers other than Z[i] and Z[a;]. 

One important difference between a general ring of imaginary quadratic integers and the Gaussian and 
Eisenstein integers is that the Gaussian integers and Eisenstein integers are not merely rings, they are 
principal ideal domains (PIDs). Hence, every ideal is generated by a singleton and one can equivalently 
work with numbers instead of ideals. The constructions of lattices over these two rings [[3l-[|5l heavily rely 
on properties of PID. However, a general ring of imaginary quadratic integers is not a PID. Therefore, in 
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order to fully exploit the potential of rings of imaginary quadratie integers, one has to work with ideals. In 
this paper, we generalize the famous Construetion A to a general ring of imaginary quadratie integers (not 
necessary a PID). We show that such construction can produce lattices that are Poltyrev-good and MSB 
quantization-good. Based on this, we show that the rates expressed in |l2|, but with function coefficients 
from the adopted ring of imaginary integers can be achieved. 

Compute-and-forward with feedback was first studied in @ in which the global channel knowledge is 
assumed at the transmitters. Using the theory of Diophantine approximation, Niesen and Whiting show 
that the traditional lattice-based scheme in [|2]| is inefficient in the asymptotically high signal-to-noise ratio 
(SNR) regime. They then proposed a novel coding scheme which is a clever combination of compute-and- 
forward and real interference alignment Q- Their scheme achieves the full degrees of freedom (DoF); 
but in order to see a gain, it requires an enormously high SNR. Another approach proposed in IS] is 
to phase-precode the lattice-based scheme in [|2l|. In this approach, one rotates the transmitted signal 
space according to the channel realization in such a way that the received signal space is close to a 
linear integer combination of the codebook. Hence, instead of the global channel knowledge, the phase¬ 
precoding approach only requires limited feedback, i.e., each transmitter only needs to know its optimal 
(or a reasonably good) phase for precoding. 

In this paper, we propose a novel framework called adaptive compute-and-forward which makes use 
of the proposed lattices. The idea is to let the transmitters adaptively choose the best ring of imaginary 
quadratic integers to work with according to the channel coefficients. It is worth noting that this approach 
only requires the knowledge of which ring the transmitters should work with and hence, limited feedback 
suffices. We show that the proposed adaptive compute-and-forward can achieve increased computation 
rates as compared to the conventional compute-and-forward scheme in Further, this can be used in 
conjunction with the phase precoding scheme in [[8]|. 

The idea of using different sets of algebraic integers for compute-and-forward was first proposed 
independently in [fT^ . In |[9l, Vazquez-Castro uses finite constellations carved from some rings of 
imaginary quadratic integers which also form Buclidean domains (hence PIDs) for compute-and-forward. 
In fTOl, instead of being confined in Euclidean domains or PIDs, we go beyond PIDs and construct lattices 
over rings of imaginary quadratic integers for compute-and-forward. However, their goodness have not 
been shown and the idea of adaptively choosing the rings of integers was only vaguely mentioned in 
idol. This paper contributes to the literature by proving the optimality of the proposed lattices, and hence, 
deriving the achievable information rates with lattices over imaginary quadratic integers. The result of this 
paper also hints at why establishing converses or capacity results for compute-and-forward problems is 
difficult as they must allow for the possibility of many rings of integers, not just Gaussian or Eisenstein 
integers. 

The main contributions of this paper are also listed as follows: 

• While most of the work in the literature using large dimensional lattices ( ifTTI - llT^ for examples) 
focus mostly on constructions over PIDs (Z, Z[i], and l^oj] in particular) and that using algebraic 
integers for lattices focus mainly on small dimensional lattices (see lfT7l - lfT9ll for example), the 
present work is one of the first in considering constructions of lattices/lattice codes over a general 
Dedekind domain, i.e., we combine the so far most popular constructions for the small and large 
dimensional lattices to construct asymptotically good lattices which inherit the algebraic structures of 
the underlying small dimensional lattices. We must emphasize that there is an independent work lf20l 
which considers lattice codes via Construction A over number fields for block fading channel. This 
construction of lattices is quite general and subsumes our construction as a special case. However, 
due to different applications at hand, they largely focus on constructions over rings of integers of 
totally real number fields and do not prove the optimality of lattices thus constructed. 

• To the best of our knowledge, this work presents the first results about the optimality of lattices 
constructed over non-PID domains. This paper advances our knowledge of lattice coding theory 
towards large dimensional random lattices generated from small dimensional ones which are full 
of structure and can serve as a stepping stone for further research. Recently, there have been some 
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researches that identify the benefits of using lattice codes constructed over rings of algebraic integers 
other than imaginary quadratic integers (see EOll and ll^ for example). The proof techniques used 
in this paper may be extendable to prove the optimality for those lattices. 

• This paper provides an application of the proposed lattices to compute-and-forward and proposes 
adaptive compute-and-forward. Simulation results show benefits of the proposed scheme which 
exploits the structure of proposed lattices and allows computation of functions with respect to the 
underlying ring of integers other than Z, Z[i], and Z[a;]. As mentioned above, Vazquez-Castro also 
proposed codes carved from lattices constructed over rings other than Z, Z[i], and Z[a;] for compute- 
and-forward in BUl. However, the rings considered therein are limited to those belonging to Euclidean 
domains, which are again PID. More importantly, no optimality has been shown in [[911 . 

A. Notations 

Throughout the paper, M and C represent the set of real numbers and complex numbers, respectively. 
We use j = \/^ to denote the imaginary unit. For a complex number x = a + jh £ £, where a, 6 G M, 
X ^ a — jb denotes its complex conjugate. We use P(E^) to denote the probability of the event E. Vectors 
are written in boldface and random variables are written in Sans Serif font. We use x to denote the 
Cartesian product and use © and 0 to denote the addition and multiplication operations, respectively, 
over a finite field where the field size can be understood from the context if it is not specified. Also, 
we do not distinguish the multiplication operation over the complex field and finite fields as it can be 
understood from the context. 

B. Organization 

The rest of the paper is organized as follows. In Section [III we discuss the compute-and-forward relay 
network which is a building block of large networks and is usually used to demonstrate the usefulness of 
the compute-and-forward technique. Construction A lattices over imaginary quadratic integers is proposed 
in Section IIII] followed by the proposed adaptive compute-and-forward scheme in Section |IVl Some 
concluding remarks are then given in Section |Vl Since the proposed construction and the proposed adaptive 
compute-and-forward heavily rely on knowledge about lattices, algebra, and algebraic number theory, a 
brief review about these topics is given in Appendix Also, we defer the proofs of the optimality of 
the proposed lattices to Appendix |Bl 


II. Problem Statement 

The network considered in this paper is the compute-and-forward relay network studied by Nazer and 
Gastpar in [|2l . The network consists of K source nodes and M destination nodes as shown in Fig [T] Each 
source node has a message Wk G {l,2,...,iy}, k G {1,... ,K} which can alternatively be expressed 
by a length-vector over some finite field, i.e., Wfc G with W = . This message is fed into 

an encoder whose output is a length-codeword x*. G C^. The codeword is subject to the average 
power constraint P per complex dimension given by 

^E||xfcf < P. (1) 

The received signal at destination m is given by 

K 

ym ^ ^ h^nk^k © (2) 

k=l 

where hmk G C is the channel coefficient between the source node k and destination node m, and 
Zm ~ CAA(0,I). One can think of this model as simply a layer in a large network; therefore, these 
destination nodes are merely intermediate relay nodes only interested in forwarding signals. Thus, instead 
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Fig. 1. A compute-and-forward relay network where Si,, Sk are source nodes and Di,..., Dm are destination nodes. 


of individual messages, eaeh destination node is only interested in recovering a function of messages 
which will be forwarded to the next layer. In the Nazer and Gastpar’s setting, functions are chosen to be 
linear combination of message^ given by 

nm = bmlQ^l® ■■■®bmKQ^K, ( 3 ) 


where bmi, • • •, bmK are elements in the same field Fp with the elements in and the operations are 
elementwise. Upon observing y^, the destination node m forms Um = Gmijm) an estimate of u^. 

Definition 1 (Computation codes). For a given equation coefficient vector = [&mi, • • •, a 

{N, N') computation code consists of a sequence of encoding/decoding functions ..., £^) / {G^, • • •, Gm) 

described above and an error probability given by 


piN) A p 


M 


IJ {u^ ^ U^} 


ym=l 


( 4 ) 


Definition 2 (Computation rate for function at relay m). For a given channel vector — [hmi, • • •, hraxY 
and equation coefficient vector b^, a computation rate R{h.m, b^, P) is achievable at relay m if for any 
£ > 0 there is an (A^, N') computation code such that 

N' > NR{h^, b^, P)/ log(p) and < a (5) 

Note that the first condition is equivalent to saying that W > 


In this paper, we consider the symmetric case where all the encoders are of the same rate. Thus, for 
a given H = [hi,...,hM] and B = [bi,...,bM], the achievable computation rate is P(H,B,P) = 
minP(hm, b^, P). Moreover, suppose there is a final destination collecting all the functions, it would be 

m 

able to recover all the messages if B is full rank. Hence, one can also define the computation rate of the 
network as follows. 


Definition 3 (Computation rate of the network). The achievable computation rate of the network is defined 
as 

fl(H,P)A WH.B.P). (6) 

B:B full rank 

Remark 4. The above definitions only consider the symmetric case in the sense that all the transmitters 
have the same power constraint P and all the encoders have the same rate. For constructions over Z 


*In general, the functions are not limited to linear combinations of messages and some other functions have been considered in 1221 
(23) for instance. 
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lattices, some asymmetric cases have been diseussed. For example, [|2l extends the eompute-and-forward 
paradigm to case where transmitters have asymmetric rates. This extension is enabled by eonstrueting 
a sequence of nested fine lattices and allowing eaeh transmitter to choose a different fine lattice in this 
sequenee. In ||2^ . Ntranos et al. further generalize the compute-and-forward paradigm to the scenario 
where transmitters may have unequal power constraints. This is done by allowing the transmitters to 
have different eoarse and fine lattiees. It is worth noting that the adaptive compute-and-forward scheme 
proposed in this paper can also be extended similarly to the asymmetric rates and asymmetrie power 
constraints case. In this paper, we only present the symmetrie result for the sake of conciseness. 

A. Open-Loop Compute-and-Forward 

In what follows, we first consider the open-loop setting where ehannel state information is only available 
at the reeeivers. In [|2]|, Nazer and Gastpar propose a novel paradigm ealled compute-and-forward where 
eaeh source node implements the same nested lattice code over Z of Erez and Zamir ffTTl and encode real 
and imaginary parts separately. This allows each relay to deeode the received signal to linear combination 
of the transmitted lattice points with eoefficients being integers and results in the following computation 
rate at a relay. 

Theorem 5 (Nazer-Gastpar). At the mth relay, given G and G the following computation 

rate is aehievableH 

a^, P) = log+ |^(^||a,^||2- ^ (7) 

Eaeh relay ean adaptively ehoose the coefficients a^ according to the ehannel vector such that 
the above computation rate is maximized. Eet A = [ai,..., a^^] and B = [bi,..., b^] be its finite field 
representative. As long as the coefficient matrix B is full rank, the eentral destination, whieh has colleeted 
all the computed functions, is able to solve the linear equations and get the individual messages. 

In [O, the real and the imaginary parts are separately considered. In what follows, we provide a 
high-level description of the scheme for the real part only but the other part works identieally. Each 
source node adopts a same nested lattice code construeted over Z. Since lattiees are closed under integer 
linear eombinations, the mth destination can now decode the codeword corresponding to an integer linear 
eombination of eodewords whieh were sent. Note that the ehannel output is a noisy version of a linear 
combination of codewords; hence, extra noise will be introduced when we try to enforee real linear 
eombinations into integer linear eombinations. This extra noise, ealled self-interference, can be equivalently 
represented as quantization error of quantizing a version of hm with respect to the quantizer Z^. After 
decoding the deeoder then maps it baek to the finite field and obtains u^. 

The proof of the above result heavily relies on two key points. The first one is the existence of good 
lattiees from this construetion. Perhaps more importantly, as recognized in (Si, the second one is that 
the mapping between Z and the finite field is a ring homomorphism so that integer combinations of 
lattice points will be eorresponding to linear eombinations over finite field. This ring homomorphism then 
allows one to map baek and forth between the finite field and real filed without ruining the strueture. In 
DU, Tunali et al. eonsidered the real and the imaginary parts jointly and generalized the compute-and- 
forward paradigm to the Eisenstein lattiees whieh are construeted from Construetion A over Z[oj\. With 
this approach, the self-interference becomes quantization error of quantizing a version of hm (now in C^) 
with respeet to the quantizer Z[uj]^. This has resulted in an inereased aehievable rate on average as Z[oj] 
approximates C better than Z[i]. 

^Note that here but in Definitionj^ computation rate is defined for hm € F^. This is not an issue by letting p the field size 

tend to oo, which is exactly what is required by the coding scheme in (2). 
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III. Lattices over Imaginary Quadratic Integers 

Since both the Gaussian integers and Eisenstein integers are rings of integers of some number fields, it 
is natural to eonsider rings of integers of other number fields. In what follows, we partieularly piek those 
rings of integers of imaginary quadratic fields. The reasons that we pick such rings are twofold. First 
of all, the ehannel coefficients we are trying to quantize lie in C, which is an extension field of M with 
degree 2. Hence, it is natural to first investigate extensions with degree 2. Seeondly, quadratic fields have 
been extensively studied and many properties have been discovered. This makes the generalization a lot 
easier. The diseussions starting from this point heavily use knowledge of lattices, algebra, and algebraic 
number theory. For baekground knowledge on these topies, the reader is referred to Appendix 

Consider an imaginary quadratie field IK = Q(\/d) with d < 0 and its ring of integer Dk = where 

f \/d, d = 2,3 mod 4, 

? - I 1 +^^ ^ ^ ^ 4 _ 1 

We now diseuss the Construetion A lattiees over Dk and some of the properties of sueh lattiees. Note 

that sinee not every Dk forms a PID, one may have to work with ideals now instead of working with 

generators of ideals as we have done in PID eases. Fet p be a prime ideal that lies above p. The norm 
of the ideal p denoted by A^(p) is equal to p^ where / G {1,2} is the inertial degree. One important 
property of Dk is that every prime ideal in Dk is also maximal. Henee, Dk/P = IFp/ and we denote by 
M. : ¥pf —)■ Dk/P a ring isomorphism indueed by this quotient ring. 

Construction A [|25ll ll26ll Fet n, N he integers sueh that n < N and let G be a generator matrix of 
an {N,n) linear eode over F^/. Construetion A over Dk consists of the following steps: 

1) Define the discrete codebook G = {x = G©y:yG where all operations are over F^/. 

2) Construct A* = A4(G) where A4 ; Fp/ —)■ Dk/P is a ring isomorphism. 

3) Tile A* to the entire to form A = A* + p^. 

Theorem 6. A is a lattice over C^. Moreover, a complex vector A belongs to A if and only if cr(A) G C 
where a = o mod p^ is a ring homomorphism. 

Proof: Since Ad is a ring isomorphism, Ad(0) = 0. Moreover, 0 G p^. Thus, 0 G A. Fet 

Ai = A4(ci) + pi, 

As = A4(C2) + P2, 

where Ci, Ci G G and pi, ps G p^. We have 

Ai + As = Ad(ci) + A4 (c2) + Pi + P 2 

= M{ci © C 2 ) + p + Pi + P 2 
= Ad(c 3 )+p 3 , ( 11 ) 

where p, p 3 G p^, C 3 = Ci © Cs G G and (a) is due to the fact that A4 is a ring isomorphism. Moreover, 
choosing Cs to be the inverse of Ci and choosing p 2 = —pi — p makes As the additive inverse of Ai. 
Therefore, A is a lattice. 

To see that A is a lattice point if and only if cr(A) G G, we note that 

A G A 

A = Ad(c) + p 
(<^^) A mod p'^ = A4(c) 

{^) (a mod p^) = c. 


(9) 

( 10 ) 


( 12 ) 
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Theorem 7. For any d < 0 square free integer, eonsider K = Q{'/d), there exists a sequenee of lattiees 
from Construetion A over Dk that is Poltyrev-good and good for quantization. 

Proof: (Sketch. Please see Appendix \B\ for details) For showing Poltyrev-goodness, we tailor the 
Minkowski-Hlawka theorem specifically for Dk and then follow the steps of Loeliger in [l27ll to show 
that with high probability, the random Construction A ensemble over Dk would produce Poltyrev-good 
lattices. For showing the MSE quantization-goodness, we modify the proof by Ordentlich and Erez 1(2^ 
where we first construct a sequence of prime ideals whose norms tend to oo for each Dk and show that 
randomly picking elements in G induces uniform distribution over (Dk/p)^- One can then follow the 
steps in [f2^ to show the MSE quantization-goodness. ■ 

Eurthermore, similar to [|29ll . one can even build multilevel lattices over Dk by Construction tta as 
follows. Eet J be a ideal whose prime ideal factorization is given by J = Fl^^p/ with p^s relatively prime. 
Erom the Chinese remainder theorem, we have 

— Ok/Pi X ... X Dk/Pl 

= F X ... X F (13) 

where fi is the inertial degree of p; in Dk- We are ready to state Construction tta over Dk- 

Construction tta Eet n’’, be integers such that if < N and let G/ be a generator matrix of an (N, n^) 
linear code over F Construction tia over Dk consists of the following steps: 

Pl 

1) Define the discrete codebook = {x = G; © y : y G (F where all operations are over F /(. 

Pi Pl 

2) Construct A* = A4(C ^,..., C^) where A4 : —)■ Dk/ 3 is a ring isomorphism. 

3) Tile A* to the entire to form A = A* + 3^. 

Similar to Construction A, one can show that Construction tt^ over Dk always produces a lattice. One 
may also show that Construction tta over Dk can produce a sequence of lattices that are Poltyrev-good 
with high probability; however, the focus of this paper is on achievable computation rates rather than 
complexity so we do not pursue this. The interested reader is referred to [[29l . 

IV. Proposed Adaptive Compute-and-Eorward 

In this section, we consider the scenario where there is limited feedback. We propose the adaptive 
compute-and-forward scheme where we use the lattices proposed in Section Un] to construct nested lattice 
codes and show that for each ring of imaginary quadratic integers, one can achieve the same rate expression 
([7]) with elements in chosen from that ring of integers. Depending on how many bits of feedback is 
available, the proposed adaptive compute-and-forward preselects a set of Dk and the idea is then simply 
to work with the ring of integers in this set that would result in the maximal computation rate. This 
approach can be used in conjunction with the phase-precoded compute-and-forward to further improve 
the performance. 

Eor a given Dk, we adapt the construction in ll28]l to construct nested lattice codes over Dk- Eet p be 
a prime ideal in Dk with N(p) = p E (2N^, 2(N^) a prime where C is a constant which guarantees the 
existence of such primes (see discussion in Appendix |B] for more details about (). Eet Ad : Fp —>■ Dk/p 
be a ring isomorphism. It is worth noting that the scheme and the result for the other class of primes 
(those stay inert) can be obtained in a similar way with slight modification of parameters. 

Eet (Cf, Cc) be a pair of nested linear code such that Cc C Cj as follows, 

Ge= {Ge0v|vGF^^}, 

G; = {G^0v|veFp™^}, 


(14) 

(15) 



where Gc is a N x rric matrix and Gf = [Gc G] with G being a N x {rrif — rric) matrix. A pair of 
(sealed) nested lattiee eodes ean be eonstrueted by the eonstruetion deseribed in Seetion |III] as 


Ac = 'lM{Cc) + 7P, 
A/ = 7 > 1 (G/) + 7 P, 


(16) 

(17) 


where 7 — y ^NP\d\ zp-i is for power constraint. We then use Aj fi V(Ac) as our nested lattice code 
whose design rate is given by 




design 


TTif — rric 

N 


log(p). 


( 18 ) 


Note that the difference between any two neighboring possible design rates is 

^ log(p) < ^ log(2CA^^) 

log(2C) + 3 log(A^) 


< 


N 


(19) 


which can be made arbitrarily small as N increases. Hence, the design rate can be set to approach any 
target rate. 

Each source node adopts a same nested lattice code Af D V(Ac) obtained by the above construction. 
Specifically, for a message G the overall encoding function acts as follows. We first pad rric 

zeros in the front to get = g(w) = and form the corresponding codeword G Cf. This 

codeword is then mapped to a fine lattice codeword t*. = 'jAA{ck) mod Ac. The transmitted signal at the 
source node k is then given by 

Xfc = (tfc - dfc) mod Ac, (20) 

where is a random dither. 

According to the channel parameters, the destination node m computes a linear combination of trans¬ 
mitted signals with coefficients = [ami, ■ ■ ■, being elements in Dk and maps this function back 

to the finite field via the ring homomorphism corresponding to the inverse of Ai. 

Specifically, the m-th receiver first forms 

K 


Ym ~ 


lYm + ^ amkdk mod Ac 


k=l 


K 


^ ^ (^mAmk^k T ^mfcdfc A Olm'^m, j mod Ac 


,fc=l 

K 


K 


^ ^ amkXk T j Of rn^m ^ ^j [^mkmk amk')^k 


mod Ac 


,fc=i 


k=l 


(i'eq,m T '^eq,m) mod Ac, 


where 


K 


f eg,m ^ ^ amk^k mod Ac, 


( 21 ) 

( 22 ) 


fc=i 


is again a lattice codeword in A/ fl V(Ac) since A/ and Ac are constructed over Dk and 


K 


'^eq,m — I 


- amfc)xfc mod Ac, 


(23) 


k=l 


is the effective noise, also called the self-interference. It then performs lattice decoding to decode to 
the nearest element in A/- to form teq,m- This estimated function is then mapped back to the finite field via 









K 




fc=l 


Fig. 2. The encoding and decoding procedure of the proposed adaptive compute-and-forward scheme. 


the canonical ring homomorphism a = M. ^ o mod p to get Um = a{teq^mh)- Note that if teq^m = ^eq,m, 
then 

fF Qmfctfc mod j = q- g^fctfc + Ag 

K 


vfc=l 


= (T I ^amfc(A^(cfc) + pfc) + A^(cc) + Pc 
\fc=i 

= ( ^ cr(amfc) O [a{M{ck)) © cr(pA:)] | © a{M{Cc)) © cr(pc 
a=i / 

^ o'(amfc) © cT (A^(cfc))^ © a{M{cc)) 

^ ^hmkQ Ck 


,k=l 


W ~ ^ ^ 


(24) 


where pfc,pc G p, Cc G Cc, (a) is due to the fact that a is a ring homomorphism, (b) follows from 
bmk = (T(amfc) G Fp, and (c) is due to the linearity of Cf and the fact that Cc G Cc © C/. Since there is 
a one-to-one mapping between the input and output of the linear code Cf, one can then invert the linear 
code Cf to get corresponding to c^. Moreover, since Cc G Cc, it only affects the first rUc positions of 
Vm- Hence, and 0^^ &m,fc © ^k lie inside the same coset. We then remove the first rUc positions via 


-,-1/ 


to get the corresponding function of messages as 

Um = g“^(Vm) = bml © Wi 


bmK © Wi^. 


(25) 


The overall encoding and decoding procedure is summarized in Fig. [2] 

Let us write = [h^i, • • •, = [omi, • • •, o-mx]^- Using the goodness results in Section Hill 

and choosing to be the MMSE estimator given by 


Q^rr) — 


1 + ||h,7^||2' 


(26) 
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Fig. 3. Comparison of computation rates for different Dk for the 2 by 2 case where /in = /i 22 = 1 and /112 = /i 2 r = y2.449. 


one can follow the steps in [|3 ll^ to show that the following computation rate is achievableEl 

Bo.(h,„.a„,P) = log+ n IKIf - j . (27) 

where log^(.) = max{0, log(.)}. i.e., as N increases, we can operate at a design rate i?design > RoKi^m, P) — 
e for any £ > 0 with vanishing pi^\ 

Note that here the subscript Dk is used to emphasize that this is obtained by working over a particular 
ring of imaginary quadratic integers Dk- Suppose we have a limited feedback of u bits, then we can pre¬ 
select A, a set of 2'^ rings of imaginary integers. To obtain the highest computation rate for the proposed 
framework, one then solves the following optimization problem to decide which Dk G ^ to work with. 

i?(H, P) = max max P£ij^(H, A, P), (28) 

invertible 

where Pok(H, A, P) = min,^ RoA^m, P)- 
A. Numerical Results 

We now provide two numerical results to demonstrate the benefits of using the proposed lattices. For 
both the results, we consider the case where there are 2 source nodes and 2 destination nodes. In Fig. [3l 
we consider fixed channel coefficients hn = 1x22 = 1 and hi 2 = h .21 = j2.449 and show the achievable 
computation rates obtained by using lattices over Dk of Q{\/d) for d G {—1, —2, —3, —5, —6, —7}. One 
can see from Fig. [3] that although Z[a;] the ring of Eisenstein integers best approximates C among all 
imaginary quadratic integers, for specific channel coefficients, it is possible that there are other rings of 
integers which have elements closer to those channel coefficients than 7j[u] does. In this example, Z[A —6] 
has elements closer to the specific channel coefficients than other rings considered in this simulation and 
hence provide the best computation rate among them in the high SNR regime. 

In Fig. m we provide average achievable computation rates by using lattices over Dk of Q(\/d) for 
d G {—1, —2, —3, —5, —6, —7}. In this figure, each channel coefficient is randomly drawn from circularly 
symmetric Gaussian distribution with variance 1. i.e., its norm has Rayleigh distribution. We average over 
10000 realizations and show that on average alternating between these 6 rings provides better performance 
than that provided by working over any of them individually. 

^We slightly abuse notation by allowing elements in am from Dk while in Definition |2] it should be bm from Fp. One can get around 
with this issue by noting that bm = cr(am) and letting p oo. 
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Fig. 4. Comparison of average computation rates for different Dk- The average is taken over 10000 pairs of channel realizations drawn 
from Rayleigh distrihution. 


B. Discussion 

We now discuss some extensions and issues for the proposed scheme. As mentioned above, one can 
incorporate the idea of phase-precoding into the proposed adaptive compute-and-forward framework to 
further improve the overall performance. Specifically, according to the channel realization, for each ring 
of imaginary quadratic integers, one can first examine the computation rate achieved by phase-precoding 
the proposed nested lattice codes constructed over that ring of integers. One can then choose the ring that 
leads to the highest computation rate to work with. Note that phase-precoded compute-and-forward with 
integers can be thought of as using rotated integers to approximate the original channel coefficients [[8l|. 
But there are algebraic integers which cannot be expressed as rotated integers; hence in general, allowing 
working over other rings of algebraic integers will result in an increased computation rate. Another way 
to exploit this is to first determine the ring of integers and then rotate the integers by phase-precoding. 
Either way will result in further reduction of the self-interference. 

A potential weakness of the proposed framework is the complexity issue. This comes from two different 
aspects. Firstly, the optimization problem in (1281) is in general very difficult to solve. Fortunately, good 
approximation algorithms have been proposed for some rings of integers (Si, JH. These algorithms are 
based on FFF lattice basis reduction algorithm [|^ . The extension of the FFF algorithms to many 
other rings of integers can be found in [l3T]| for example. Moreover, simulation results shown above 
suggest that one does not have to consider too many Dk for getting improved performance. Secondly, the 
modulo operation with respect to an ideal may cause increased complexity. Fuckily, there are algorithms 
available which have polynomial running time. For example, Algorithm 1.4.12] will produce a unique 
canonical coset representative very efficiently. Furthermore, for those Dk which also form Euclidean 
domains (there are exactly five of them corresponding to d G {—1, —2, —3, —7, —11}), one can further 
reduce the complexity by taking advantage of Euclidean functions as reported in [0 . 

V. Conclusion 

In this paper, we have moved beyond PIDs and generalized Construction A of lattices to general rings 
of algebraic integers of imaginary quadratic fields. We have then shown that such construction can produce 
good lattices in the sense of Poltyrev and MSE quantization. When used for compute-and-forward, these 
lattices have allowed us to reliably compute linear combinations of codewords with coefficients being 
elements in the underlying ring which the lattices are constructed over. A novel scheme named adaptive 
compute-and-forward has been proposed where the transmitter first chooses a ring of algebraic integers 
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depending on the ehannel eoeffieients and then uses a lattiee eode over the ehosen ring of algebraie 
integers. This allows us to obtain higher eomputation rates than using a fixed lattiee eode. Moreover, one 
ean phase-preeode the proposed adaptive eompute-and-forward seheme to further improve the performanee. 

Appendix A 
Preliminaries 

In this appendix, we review the literature on lattiees, algebra, and algebraie number theory that will 
be the foundation of this work. All the Lemmas are provided without proofs for the sake of brevity; 
however, their proofs ean be found in standard textbooks. For details, please see for example [fTTI . [|T^ . 
m, 1133-1351. 


A. Lattices 

An A^-dimensional lattice A is a discrete subgroup of which satisfies the following properties: 
0 G A, VA G A, we have —A G A, and VAi, A 2 G A, we have Ai + A 2 G A. Some important operations 
and notions for lattices are defined as follows. For a x G the nearest neighbor quantizer associated 
with A is denoted as 

Qa(x) = a G A; ||x — A|| < ||x — A'll VA' G A, (29) 

where ||.|| represents the L 2 -norm operation and the ties are broken arbitrarily. After defining lattice 
quantization, we can define the fundamental Voronoi region Va as Va = {x : Qa(x) = 0}. The mod A 
operation simply provides quantization error with respect to A and is represented as 

X mod A = X — Qa(x). (30) 


The second moment of a lattice is defined as 


a\K) 


1 1 
VoI(Va) ^ 



(31) 


where VoI(Va) is the volume of Va and the normalized second moment of the lattice is then defined as 


G(A) 


Vo1(Va)2/^' 


(32) 


The normalized second moment is a dimensionless quantity and is invariant to scaling. A lower bounded 
on G(A) can be obtained by the normalized second moment of a ball which asymptotically approaches 
in the limit as N -^ 00 . 

ZTre 

We are ready to define the two goodness properties considered in this work, namely the goodness for 
MSB quantization and goodness for channel coding (Poltyrev-goodness). 


Definition 8 (Goodness for MSB Quantization). A sequence of lattices is asymptotically good for MSB 
quantization if 

hm G(A) = (33) 

n^oq ine 


Consider the communication channel y = x + z where each element in z is i.i.d.~ A/^(0, rf) and there 
is no power constraint on x G A. The Poltyrev-goodness is defined as follows. 


Definition 9 (Poltyrev-Goodness). A sequence of lattices is asymptotically Poltyrev-good if whenever 


, V°I(Va) 

27re ’ 


( 34 ) 


the error probability of decoding x from y can be made arbitrarily small. 
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B. Algebra 

Let 7?. be a commutative ring. Let a, 6 7 ^ 0 G 7?. but a 6 = 0, then a and b are zero divisors. If 
ab = ba = 1, then we say a is a unit. Two elements a,b G TZ are associates if a can be written as the 
multiplication of a unit and b. A non-unit element r G 7^ is a prime if whenever r divides ab for some 
a, 6 G 7^, either r divides a or r divides b. An integral domain is a commutative ring with identity and 
no zero divisors. An additive subgroup 3 of TZ satisfying ar G J for a G J and r G 7?. is called an 

ideal of TZ. An ideal IT of 7^ is proper if 3 ^TZ. An ideal generated by a singleton is called a principal 

ideal. A principal ideal domain (PID) is an integral domain in which every ideal is principal. Famous 
and important examples of PID include Z, Z[i] and Z[u]. Let a,b &7Z and 3 be an ideal of TZ', then a is 
congruent to b modulo 3 if a — b E 3. The quotient ring 7^/3 of 7?. by is the ring with addition and 
multiplication defined as 

{a + 3) + {b + 3) = {a + b) + 3, and (35) 

{a + 3) ■ {b + 3) = {a ■ b) + 3. (36) 

A proper ideal p of 7?. is said to be a prime ideal if for a,b E TZ and ab E p, then either a G p 
or 6 G p. A proper ideal 3 of 7?. is said to be a maximal ideal if 3 is not contained in any strictly 
larger proper ideal. It should be noted that every maximal ideal is also a prime ideal but the reverse 
may not be true. Let TZi, TZ 2 , ..., TZl be a family of rings, the direct product of these rings, denoted by 
7^1 X 7^2 X • • • X T^L.is the direct product of the additive abelian groups TZi equipped with multiplication 
defined by the componentwise multiplication. 

Let 3i and 32 be two ideals of TZ, we shall now define some operations of ideals. The sum of two 
ideals is the ideal defined as 

Ji J 2 — {ti b : a E Ji,b E J 2 } • (37) 


Two ideals are relatively prime if 7?. = 3i + 32- The product of two ideals is the ideal defined as 


rv ^ A 
^ 1^2 — 


aibi : ai E 3i, h G 32, n G 


N 


i=l 


(38) 


In general, we have 3i32 C 3i fi 32; but if they are relatively prime, then 3i32 = 3i fi 32. We say 3i 
divides 32 or 3i|32 if there is an ideal 33 such that 32 = 3i33 (this is equivalent to 32 C 3i). 

Let TZl and TZ 2 be rings. A function a : 7^i ^ 7^2 is a ring homomorphism if 


a{a + b) = a(a) 0 a(b) Va, b eTZi and (39) 

a{a ■ b) = a{a) 0 crib), Va, b eTZi, (40) 

where (0, •) and (0,0) are operations in TZi and TZ 2 , respectively. A homomorphism is said to be 
monomorphism if it is injective and isomorphism if it is bijective. Let TZ he a commutative ring, and 
3i,..., 3„ be ideals in TZ. Then, from the Chinese Remainder Theorem, we have 


TZI 3i = 7Z/3i X ... X 7^/3„. 


(41) 


C. Algebraic Numbers and Algebraic Integers 

Now, we provide some background knowledge on algebraic number theory. The materials are mostly 
from [l34l [f35ll and proofs and algorithms can be found therein. 

Definition 10 (Algebraic Numbers and Algebraic Number Fields). An algebraic number is a root of some 
polynomial with coefficients in Z. The set of all algebraic numbers is a subfield A of C. We define a 
number field to be a subfield IK of A (hence a subfield of C) such that the degree [K : Q] is finite. 

Theorem 2.2 in If34l shows that any such IK is equal to 0(0), the smallest subfield containing Q and 
9, for some algebraic number 9. 
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Definition 11 (Algebraic Integers). An algebraic integer is a complex number which is a root of some 
monic polynomial (whose leading coefficient is 1) with coefficients in Z. The set of all algebraic integers 
forms a subring B of C. For any number field K, we write Dk = and call Dk the ring of integers 

of K. 

From Corollary 2.12 in [|34l . one has that if K is a number field then K = Q{9) for an algebraic 
integer 9 which is called a primitive element for K over Q. Also, in general, there will be several distinct 
Q-monomorphisms (i.e., it fixes Q) embedding IK into C. From Theorem 2.4 in [[341, we know that for 
IK = Q{9) a number field of degree n over Q, there are exactly n distinct Q-monomorphism cr* : IK — )■ C 
and such monomorphisms form a group Gal(IK/Q) = {ai, ... ,an} which is referred to as the Galois 
group. Moreover, for a G Q{9), ai{a) for i G {1,2,... ,n} are the distinct zeros in C of the minimal 
polynomial of a over Q. We call those crj(a) the conjugates of a and define the norm of a to be the 
product of conjugates as 

n 

^K(a) = JJo-i(a)- (42) 

i=l 

Let {«!,..., an} be a Q-basis for K. We define the discriminant of {ai,..., an} as 



/ ai(Q!i) 

cri(a2) . 


\ 

A[ai,..., an] — det 

(T2(Q!i) 

0'2(q<2) • 

• 0'2{<^n) 



V <7n(Q;i) 

<7n(Q<2) • 


/ 


If {«!,..., an} is a Z-basis for Dk, we define the discriminant of K to be Ak = A[ai,..., an]- Let 
3 be an ideal o f Dk- The norm of 3 is N{3) = IDk/UI- Moreover, if {^i,... ,/9n} is a Z-basis for 3, 

then N{3) = The norm is multiplicative, i.e., for two ideals 3i and 32 of Dk, N{3i32) = 

N{3i)N{32). It can be shown that if N{p) is a rational prime, then p is a prime ideal. 

In this paper, we are particularly interested in imaginary quadratic fields and their algebraic integers 
whose definitions can be found below. 


Definition 12 (Quadratic Fields). A quadratic field is an algebraic number field IK of degree [IK : Q] = 2 
over Q. Particularly, one may write K = Q{Vd) where d G Z is square free. We say IK is an imaginary 
quadratic field if d < 0. 

Let K = Q(\/d). One can show that Dk = '^[i] where 

J \/d, d = 2, 3 mod 4, 

d^l mod4. 

Also, Ak = 4d if d = 2, 3 mod 4 and Ak = d if d = 1 mod 4. It can be easily seen that when d = —1 
we have the Gaussian integers and when d = —3 we have the Eisenstein integers. 

Example 13. Let us consider the case d = —5, i.e., K = Q(\/^). Let a = a + where a,b e Z. 

Since the degree is 2, there are exactly two Q-monomorphisms. In order to have a Q-monomorphism, 
one must have a{^/—5) ■ cr(\/^) = cr(—5) = —5, which implies that = ±\/—5. Thus one has 

that 

ai{a) = a + bai{-\/—b) = a + by/—5, 

<72 (tt) = a + ba2{y/—5) = a — by/—5. (45) 

Then the norm of a is ai{a) ■ a 2 {a) = + 56^ which coincides with the Euclidean norm. Since —5 = 3 

mod 4, from ([8]), {1, y/^} is a Z-basis for Dk- One can calculate the discriminant as follows, 

Ak = det J ^ = -20. (46) 
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For any prime p, pX is a prime ideal in Z. Let p be a prime ideal of Dr. We say p lies above p if p \pX. 
The ideal pO^ can be uniquely faetorized into with p/ distinet. We eall ei the ramification 

index of p; over p and fi = [Dk/Pz : Z/pZ] the inertial degree of p; over p. Note that one must have 
N{Vi) = Also, the ramifieation indiees and inertial degrees must satisfy = n. If ei > 1 for 

some I, we say p (or pD^ to be precise) ramifies in Dk- If L > 1, we say p splits in Dk- If L = 1 
and Cl = 1 (i.e., /i = n), we say p remains inert in Dk- The following Lemma allows one to efficiently 
categorize primes. 

Lemma 14. Let p be a rational prime. For a quadratic field K, one has 

• if fj = 0, then p ramifies in Dk, 

• ^ P split® 

• if f = — 1, then p remains inert in Dk, 

where is the Kronecker symbol mod p. Moreover, the Kronecker symbol mod p operation can 

be efficiently computed. (See for example |l36l Algorithm 1.4.10].) 

It is well-known that not every ring of imaginary quadratic integers forms a PID. In fact, it was 
conjectured by Gauss and shown by Heegner and Stark that there are only 9 of them that are PIDs 

(corresponding to d G { — 1, —2, —3, —7, —11, —19, —43, —67, —163}). Therefore, it is crucial to have a 

systematic way to identify prime ideals. The following theorem provides a means of doing this. 

Lemma 15. Let p be an odd rational prime. For a quadratic field K = Q(\/d), one has 
. if p ramifies in Dk, then p = (p, y/d) is a prime ideal lying above p, 

• if p splits in Dk, then p = (p, a + y/d) is a prime ideal lying above p for any a such that = d 

mod p. 

Moreover, such a can be efficiently found (See for example [1^ Algorithm 1.5.1].) 

One important property of Dk is that every prime ideal is maximal. Therefore, we have that for every 
prime ideal p in Dk, 

1^k/P —IFp/, (47) 

where / is the inertial degree described above. 

Example 16. Again consider d = —5, i.e., IK = Q(\/^), and p = 23. We have that = 1; hence, 

23Z splits into two prime ideals in Dk- From the above theorem, one can check that 8^ = —5 mod 23; 
thus, Dk = PP where p = (23,8 + y/^). Also, we have 

Therefore, A^(p) = y/—10580/ — 20 = 23. Moreover, Dk/P — IF 23 - This coset decomposition and the 
corresponding ring isomorphism is shown in Fig. [S] 

Lemma 17 (Dirichlet’s prime theorem). For any two relatively prime integers a and d, there are infinitely 
many rational primes of the form p = a mod d. 

Note that for any odd prime p, depends on p mod 4Ak. Hence, Lemma [HI and Lemma [TtI 

together imply that for any quadratic field K, there exist infinitely many splitting primes and infinitely 
many inert primes as well. Moreover, for Dk of Q{y/d) with d = 2, 3 mod 4, the behavior of its primes 
is regulated by 
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Fig. 5. The coset decomposition and the corresponding ring isomorphism. 


where we have used the multiplicative property of Kronecker symbol and the fact that f ^ j = 1. Therefore, 
for this class of rings, the behavior of the prime p only depends on p mod dd. In what follows, we provide 
a weaker version of the Chebotarev’s density theorem which further tells us how those primes distribute 
asymptotically. 

Lemma 18 (Chebotarev’s density theorem). In a ring of algebraic integers of a quadratic field, asymptot¬ 
ically, the density of each category of primes is 1/2, i.e., asymptotically, half of the rational primes split 
and half of them remain inert. 


Appendix B 
Proofs 

In this appendix, we show that Construction A over Dk produces good lattices with high probability. 
Throughout the proof, we only consider Dk, the ring of algebraic integers for Q(\/d) with d = 2, 3 
mod 4; the case of d = 1 mod 4 can be proved similarly with a slight modification of parameters. Also, 
we focus on those primes which split completely in Dk for the sake of simplicity; however, the proof 
techniques can be applied to the proposed construction with primes that remain inert. 


A. Poltyrev-Goodness 

The proof closely follows the steps in llTTll . Let p be a splitting prime in Dk, i.e., pDk splits into two 
prime ideals in Dk, namely pD^ = pp. Therefore, we have Dk/p = Dk/P = Fp. Let C be the collection 
of all {N, n) linear codes C over Fp. The set C is a balanced set and the basic averaging lemma llTTl 
Lemma 1] applies. Thus, one has 


' ' ceCcec\o 


p" - 1 

pX _ 1 


se(F^)\o 


(50) 


for an arbitrary mapping / : —?■ M. Since we can identify C by the basic averaging lemma works 

for arbitrary mapping / : —)■ M as well and we use and interchangeably in the following. 

Theorem 19 (Modified Minkowski-Hlawka Theorem). Let / : M be a Riemann integrable function 

of bounded support. Then, for any integer 0 < n < A^, and any fixed VoI(Va), the approximation 


1 


~ voi(v.yA) ^ 

CgC vG7A\0 


/(v)dv, 


(51) 
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becomes exact in the limit p —)■ oo, ^ 2 ^ vI^kI ^ VoI(V^a) = "■ fixed. 

Before proceeding to the proof, we first note that due to Dirichlet’s prime theorem (Lemma [T7]) and 
Chebotarev’s density theorem (Lemma [TSl) . there exist infinitely many splitting primes in every Dk so 
that one can safely let p go to infinity. 

Proof: Recall that a = o mod p^. Note that by the basic averaging lemma, 


' ' ceCve 7 A\o ' ' cec 




ve(D^)\0:(T(v)=0 


ve(OK 


/(7v) 

v6(i)^)\0:cr(v)=0 

/(7v) 

ve(O^)\0:<7(v)=0 


EE 

CeC ceC\o 


/(7v) 


veO^ :tT(v)=c 


p^-l 

+ - 1 2-^ 


pN _ I 


Y1 + 

ve(i)^)\0:(T(v)=0 






ceF^\o 
p" - 1 
jgiV _ 1 

-N 


Y1 


veO^ :o-(v)=c 


/(7V) 


veO^ :cr(v)7^0 


(b) 


VoI(V^a) 


ve(7OK)^:o-(v)^0 


E /('') 7 '' { 




N 


(52) 


where (a) requires 7 p being large and / having bounded support and (b) requires 7 ^( ) to be small 

so that the Riemann sum approaches the Riemann integral. ■ 

One can then follow the proof in [[27l to show that with high probability, the proposed construction 
produces lattices that are Poltyrev-good. 


B. MSB Quantization-Goodness 

The proof closely follows the steps in lf2^ . We again consider splitting primes that can be decomposed 
into p = pp and use the proposed construction with the prime ideal p. Denote by Vat the volume of an 
A^-real dimensional ball with unit radius and let i3(s, r) be a 2A^-real dimensional ball in (or equivalently 
A^-complex dimensional ball) with radius r centered at s. We again prove the result for p primes splitting 
completely in Dk only, i.e., A^(p) = p for a prime ideal p lying above p. Note that scaling would not 
change lattice structure; in the sequel, we equivalently consider the scaled version 

K = ^M{C)+^p^, 


(53) 
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where 7 = \jANP\d\ zp-i and we piek 2iV^ < p < 2(^N^ where C is a eonstant to guarantee that we 
ean find a prime p splitting in Dk- Note that as mentioned in Appendix the behavior of primes in 
imaginary quadratic integers is regulated by the Kronecker symbol j and depends only on p mod id 
for the class of rings considered here, i.e., we can find at least one a such that every natural prime p 
in the arithmetic progression p = a mod D splits completely. The prime number theorem for arithmetic 
progressions llTTl states that for any > 1 , there must exist a sufficiently large c such that for every 
z > c, there’s at least one p G {z, (z) in this arithmetic progression. However, it may not be an easy 
task to characterize how large c has to be. Fortunately, the results in IfJTl Table I] provide small upper 
bound on c with small ( for many cases of D < 840, which is far more than enough for the application 
of adaptive compute-and-forward. 

Lemma 20 (Modified Lemma 1 in |[28l). For any s G and r > 0, the number of points of inside 
B{s,r) can be bounded as 


niax{r - 0} 


2N 


V 2 N 


< l-^K nS(s,r) 


< r + 




v I^kI 

2 

27V 


N 


V 2 N 




TV ■ 


(54) 


Proof: Let S = fl r) + V£, N where is the funci 3 .ment 3.1 Voronoi region of One h 3 .s 

K K 


Vol(5) = |£)Jne(s,r)| 



(55) 


Note that the largest distance between any two points lying within Vo^v can be upper-bounded by ^j2N\d\. 
For any x G i3(s,r — ^J2N\d\/2) that lies within a + V^iv for some a G we have 


la — x|| < 


y^2N\d\ 


(56) 


Therefore, following the triangle inequality, we have ||a — s|| < r and thus x G 5. This shows x G 
B{s,r — ^j2N\d\/2) C S. Moreover, we have 


5 C S(s, r) + Vqn C B(s, r) + S(0, ^27^/2), 


(57) 


which shows that S C B{s,r + ^J2N\d\/2). Evaluating the volumes of these sets completes the proof. 
For any x G C^, define 

dfx, A) =-min ||x — A||^ 

^ ’ 2N ASA " " 


=- min ||x — 7 A 4 (c) — 7 a|P 

2Araep^,cec" ' ^ ^ " 


-min ||(x — 7 A 4 (c))*|P, 

2N cec ' V 77 II > 


(58) 


where y* = y mod 7 p^. Also, note that 


d(x, A) < 


^ 7 


< 




27r 


(59) 
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where is the minimum squared Euelidean distanee of elements in p and the lower bound is due from 
a bound on dp 2 [[38l . [[391. Reeall that for the ease eonsidered (d = 2,3 mod 4), Ak = 4d. For any 
w G Fp \ 0, define the random veetor C(w) = G © whieh is uniformly distributed over F^. Thus, 
7 V1 (G(w)) is uniformly distributed over (Dk/p)^. 

For all w G Fp \ 0 and x G C^, we have 


£=P(^ll(x-7X(C(w)))-f <:f 


= p 
= p 

(a) 

> p 


-N 


-N 


7(0^^)"^ ns* (x, Vnp) 


7l)^nS(x, VNP) 


-N I ^-1 


^VNF- 


V^W\ 


2N 


V 2 N 




N 


= V2N{l-''NP\d\-^p-^)^ I 1 


7 


VWl ] 

J 


2N 


A) T/ o-2Ar 
— V2N2^ 


(c) 

> V2n2-^^ 


1 - 


^y2N\d\y^' 


2N 


Vp 


N 


2N 




{2Ny 


( 60 ) 


where (a) is from Femma[20l (b) is due to the ehoiee 7 = y 4AP|d| 2 , (c) is due to the ehoiee p > 2N^, 
and (d) is true for suffieiently large N. This ean be verified by noting that (1 — \d\^/^/N)'^^ is a positive 
strietly inereasing funetion for N > \d\^/^ and will eonverge to exp(—2|(i|^/^) > 0 while 1/(2A)^ is a 
positive strietly deereasing funetion and will eonverge to 0. 

Fet W = p^ — 1 and label eaeh of the wGFp\0byi = Define the indieator random 

variable related to x G as 

./l, ^\\(yi-rp-''^M(Cm?<P 

0 , otherwise 


Xi = 


( 61 ) 


One has that for any x G 

P((d(x,A))>P) = P( 5 ^X. = 0 


w 


i=l 

W 


< P 


i=l 
1 


> £ 


, w w 




EE Cov(xi, Xi) 


i=l 1=1 


< ^ < C(2Nfp-’^2^^VfZ, 


(62) 























20 


where (a) follows from the Chebyshev’s inequality, (b) is due to the faet that C(wi) and C'(w 2 ) are 
statistieally independent unless wi = a ■ W 2 for a G Fp, and (e) is by plugging in (1^ and the faet that 

P < 2(N^. 

One ean then show that for any distribution on X, we have 


Ex,A(d(X,A)) 



1^1 + 5{2Nf2 




where 5 = ('/tt is a eonstant. This in turn implies that 


Jim Ea (cr^(A)) < — per real dimension, 
if one ehooses the eoding rate to be 


n 

N 


log(p) = log 



+ e, 


(63) 


(64) 


(65) 


for e > 0. The volume of the fundamental Voronoi region is lower bounded by 


Henee, we have 


/ N \ 2/2iV 

= 2-^NPV^l^^^. 


lim Ea (Ga) = lim Ea ( 

2N^oo 2N^oo \Vo1(Va)^/^" / 

EA(a 2 (A)) 


^ 2iV™oo 2-^2Arp^2j2iv 

= 2" lim 

2N^oo2N 


= T 


1 


27 rexp(l) 


( 66 ) 


(67) 


Similar to [|28l . one ean then use the above result to show that asymptotieally, most of the lattiees thus 
construeted will be good for MSE quantization. 
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